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ABSTRACT 

We show analytically that large abundances of NH2D and NHD2 can be produced by gas phase 
chemistry in the interiors of cold dense clouds. The calculated fractionation ratios are in good agreement 
with the values that have been previously determined in L134N and suggest that triply-deuterated 
ammonia could be detectable in dark clouds. Grain surface reactions may lead to similar NH2D and 
NHD2 enhancements but, we argue, are unlikely to contribute to the deuteration observed in L134N. 
Subject headings: ISM: individual objects: L134N - ISM: chemistry - ISM: deuterium 



1. INTRODUCTION 

Recently NHD 2 and NH 2 D were detected in L134N 
by Roueff et al. (2000). This is the first detection of a 
doubly-deuterated molecule in a dark interstellar cloud. 
Roueff et al. derived fractionation ratios of i?(NH2D) = 0.1 
and i?,(NHD2) = 0.05, where we define the fractiona- 
tion ratio, R, of the deuterated species XH m Dm' to be 
n,(XH m D m /)/n(XH m+ iD m '_i) where n is the number den- 
sity. NH2D was previously detected at the same position 
by Olberg et al. (1985) and Saito et al. (2000), who deter- 
mined i?(NH2D) ss 0.05. The first interstellar detection of 
N2D + was also made at approximately the same location 
by Snyder et al. (1977), who derived R(N 2 D+) = 0.45. 
DCO + has been detected in L134N by Wootten, Loren, & 
Snell (1982), Guelin, Langer, & Wilson (1982), and But- 
ner, Lada, & Loren (1995) with values for i?(DCO + ) in 
the range 0.03-0.07. More recent observations of N2D" 1 " 
and DCO + at the same location as the NHD2 peak were 
carried out by Tine et al. (2000), who obtained fractiona- 
tion ratios of 0.35 and 0.18 respectively for these ions. In 
addition to the deuterium fractionation, it is known that 
the absolute abundances relative to H2 of NH3 and N2H 4 " 
peak at the same position (Ungerechts, Walmsley, & Win- 
newisser 1980; Swade 1989; Dickens et al. 2000). 

The only other doubly deuterated molecule observed 
in the interstellar medium is D2CO, which has been de- 
tected in two regions: the well-known Orion Compact 
Ridge source (Turner 1990) and the low-mass protostar 
IRAS 16293-2422 (Ceccarelli et al. 1998; Loinard et al. 
2000). In both cases, the large D2CO fractionation - 
i?,(D2CO) = 0.02 and 0.35 respectively - is thought to 
result from grain surface chemistry, since the molecular 
abundances in these warm regions reflect the evaporation 
of ice mantles from interstellar dust grains (e.g. Brown, 
Charnley, & Millar 1988), and the post-evaporation molec- 
ular D/H ratios remain equal to the ratios in the precursor 
ices for over 10 4 years (Rodgers & Millar 1996). Deuter- 
ated formaldehyde forms on grains via addition of H and D 
atoms to CO (Tielens 1983; Charnley, Tielens, & Rodgers 
1997), and so large amounts of HDCO and D2CO can be 
expected if the gas-phase atomic D/H ratio is large. 

NH2D and NHD2 can also form on grains through D 
and H atom additions to atomic N (Brown & Millar 1989). 
However, as noted by Roueff et al. (2000), the lack of on- 



going star formation, together with the low temperature 
in L134N (T w 9-12 K; Swade 1989; Dickens et al. 2000), 
suggests that mantle removal has not occurred. Markwick, 
Millar, & Charnley (2000) have proposed a nonthermal 
mechanism for removing grain mantles in dark clouds in- 
volving ion-neutral streaming produced in MHD motions. 
However, this process also acts to remove H2D + and N2D" 1 " 
from the gas (Charnley 1998) and so this mechanism is un- 
likely to have occurred at the deuterium emission peak in 
L134N. It therefore appears that the NH 2 D and NHD 2 
observed in L134N must be created in the gas phase. 

In this paper, we address the issue of gas phase ammo- 
nia deuteration. We first discuss the underlying chemistry 
which controls the deuterium fractionation and ammonia 
abundances in dark clouds. We then show how successive 
deuteron transfer reactions can lead to large abundances 
of multiply-deuterated ammonia. We derive analytical ex- 
pressions for the steady-state D/H ratios in isotopomers of 
ammonia, and compare these theoretical ratios with those 
observed in L134N and with those expected from grain 
surface chemistry. 



2. DEUTERIUM FRACTIONATION IN DARK CLOUDS 

A number of singly-deuterated molecules have been de- 
tected in dark interstellar clouds with large fractionations 
(R ~ 0.001-0.1), and the theory behind the observed D en- 
hancement is well understood (e.g. Watson 1974; Guelin 
et al. 1982; Millar, Bennett, & Herbst 1989; Roberts & 
Millar 2000). Essentially, small zero-point energy differ- 
ences ensure that molecular ions become preferentially 
deuterated via D/H exchange reactions with HD, and sub- 
sequent ion-molecule reactions spread this D-cnrichment 
to neutral species. At 10 K, the most important fraction- 
ation process is that of H2D+: 
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We can quantify the degree of fractionation attainable by 
reaction (fil) by writing 



i?(H 2 D+) = SR(HD) 



(2) 



where i?(HD) will be twice the cosmic D/H ratio of ~ 
1.6 x 10~ 5 (Linsky et al. 1995), and <S is the so-called en- 
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hancement factor, given by: 



S 



^2jkjx{rrij) + a e (H 2 D+)a 



(3) 



(Stark, van der Tak, & van Dishoeck 1999; Millar et al. 
2000). kf and k r are the forward and reverse rate coef- 
ficients for reaction ([I]); at 10 K the respective values are 
1.7 x 10~ 9 and 3.6 x 10" 18 em's" 1 (Millar et al. 1989; 
Sidhu, Miller, & Tennyson 1992), x e is the electron frac- 
tion, kj is the rate coefficient for proton or deuteron trans- 
fer to species rrij (principally CO, N 2 and O) whose abun- 
dance is x(m,j) = n(m./)/n(H 2 ), and a e (H2D + ) is the elec- 
tron dissociative recombination coefficient of H2D + (equal 
to 6 x 10~ 8 (T/300K)- ' 65 cm 3 s -1 ; Larsson et al. 1996). 

Proton transfer from to CO and N 2 is the primary 
source of the HCO + and N2H + ions, so the H2D + en- 
hancement also determines the fractionation of DCO + and 
N2D+. If we assume that one third of such reactions result 
in deuteron transfer, it follows that 

i?(DCO+) = i?(N 2 D+) = ±i?(H 2 D+) = §Si2(HD) (4) 

This simple theory predicts that the fractionation of N 2 D + 
and DCO + should be equal. However, as discussed in §1, 
the values measured in L134N by Tine et al. (2000) are 
0.35 and 0.18 respectively. This discrepancy is puzzling, 
since if anything, the value of i?(DCO + ) should be larger 
than i?(N2D + ) as DCO + is also fractionated via the reac- 
tion of HCO+ with atomic D (Adams & Smith 1985; cf. 
Dalgarno & Lepp 1984; Opendak 1993). A possible expla- 
nation is that deuteron transfer from H2D" 1 " to N 2 occurs 
preferentially to proton transfer, but that the branching 
ratios for H 2 D + + CO are statistical. Nevertheless, the 
N2D+ fractionation is comparable to that determined by 
Snyder et al. (1977) and, to order of magnitude, the ob- 
served fractionation is compatible with the above theory 
if S > 10 4 . 



at 10 K, chemical cycling between ortho and para-H 2 via 
proton transfer reactions is able to provide sufficient ortho- 
H 2 (Le Bourlot 1991, 2000). 

At the densities appropriate for L134N, gas phase chem- 
istry, rather than freeze-out onto grains, is the dominant 
loss route for ammonia. Cosmic ray ionization of H2 fol- 
lowed by successive proton transfer reactions results in 
NH3 molecules being recycled into NH|" ions, which can 
then recombine to either NH3 or NH2. In the latter case, 
the NH2 is rapidly destroyed by atomic oxygen. The 
resulting expression for the steady-state ammonia abun- 
dance is: 



z(NH 3 ) = 



0.3eCa;(Hc)a;(N2) 
(l-eMCO) + x(N 2 )]r H 4 



(7) 



where e is the fraction of NH4 dissociative recombinations 
that lead to NH3, ( is the cosmic ray ionization rate, and 
T H + is the total rate (s _1 ) for proton transfer to NH3. In 
theory, T H + will depend on the abundances of all molecu- 
lar ions capable of transferring a proton to ammonia, but 
because HCO + , N2H + and H^ are the dominant ions we 
can write 



10~ 9 [2.2n(HCO+) + 2.3n(N 2 H+) + 2.7n(H+)] 



(8) 



Based on the fact that the abundances of HCO + and N2H" 1 " 
are observationally well-determined in L134N, and that 
the abundance of H.j~ can be obtained from simple chem- 
ical arguments (see §H), we derive T H + pa 10~ 12 s _1 . The 
value of e was measured in the laboratory by Vikor et al. 
(1999), and was found to be 0.69. With assumed values 
for a;(He) = 0.15, and ( = 1.3 x 10~ 17 s -1 , the observed 
ammonia abundance in L134N can be reproduced by equa- 
tion (Q) if the N2/CO ratio is ps 0.1, consistent with the 
N2 abundances determined in a number of dark clouds by 
Womack, Ziurys, & Wyckoff (1992b). 



3. GAS PHASE AMMONIA FORMATION 

The observed ammonia abundance in L134N is ~ 10~ 7 
(Swade 1989; Dickens et al. 2000), in agreement with the 
steady-state value in chemical models where most of the 
nitrogen is in molecular form (e.g. Millar, Farquhar, & 
Willacy 1997). In this case, ammonia formation is initi- 
ated by He + attack on N2 to form N + . This ion then un- 
dergoes successive reactions with H2 until NH|" is formed, 
which then recombines to give NH3: 



N 2 ^ NT 



4H, 



NEC 



NH, 



(5) 



This formation mechanism remains somewhat controver- 
sial, because the first hydrogenation step 



N 



H 2 



NH 



H 



(6) 



is endothermic by ps 170 K (Marquette, Rebrion, & Rowe 
1988). However, this is slightly less than the ground 
state rotational energy of ortho-H2, and Le Bourlot (1991) 
showed that as long as the ortho-to-para ratio is greater 
than ~ 10~ 4 then reaction (|q) is the dominant loss route 
for N + , so reaction sequence (||) proceeds efficiently. Even 



4. GAS PHASE AMMONIA FRACTIONATION 

We have extended our steady-state analysis to include 
multiply deuterated ammonia. Figure 1 shows the pri- 
mary chemical reactions creating and destroying NH3 and 
its associated isotopomers through ion-molecule reactions 
involving generic protonated and deuterated ions, XH + 
and XD + . After NH3 is formed by the reaction sequence 
(|^), deuteron transfer reactions form NH3D + which can 
then recombine to give NH 2 D. Successive deuteron trans- 
fer reactions can lead eventually to NHD2 and ND3. The 
relative steady-state abundances therefore depend on the 
XD+/XH+ ratio and the branching ratios for dissociative 
recombination of the deuterated ions. 

We can quantify the XD + /XH + ratio by introducing the 
parameter R, equal to the ratio of the rates for deuteron 
vs. proton transfer, i.e. we define R = r D +/r H +, where 
r H + is given by equation (|^) and r D + is obtained from 

T D+ re 10~ 9 [2.2?i(DCO + ) + 2.37i(N,D+)+0.9n(H 2 D + )] 

As with r H + , there will also be a small contribution to 
T D + from isotopomers of less abundant ions (e.g. H3O" 1 ", 
HOCO+, HCNH+, etc.), but because HCO+ is the most 
abundant ion, and because the deuterium fractionation in 
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Fig. l. — Reaction network for ammonia deuteration in cold clouds. XH + and XD + represent all species capable of 
transferring a proton or deuteron to NH 3 ; principally , N 2 H + , HCO + , and their deuterated isotopomers. Also shown 
are the assumed statistical branching ratios for each dissociative recombination channel, e is the fraction of dissociative 
recombinations which result in an isotopomer of NH3 rather than NH2 . 



these minor ions originates with deuteron transfer from 
H 2 D+ and DCO+, in practice R w i?(DCO+). 

We assume that the NH3/NH2 branching ratio, e, is the 
same for deuterated ammonium ions as for NH^, and we 
assume statistical branching ratios regarding the position 
of the D in the products (see Fig. 1). Gellene & Porter 
(1984) measured the branching ratios for electron dissocia- 
tive recombination of NHDjj" and found H atom ejection 
occurred 2.7 times more frequently than D ejection, im- 
plying that the N-H bonds are eight times more likely to 
break than the N-D bond. A similar effect has been ob- 
served in the recombination of HDO + and HD20 + (Jensen 
et al. 1999, 2000). Therefore, our results represent a lower 
limit to the amount of deuteration that can occur via ion- 
molecule chemistry, since preferential retention of the D in 
the molecule is more likely to occur. 

For each of the species in Fig. 1, one can write an expres- 
sion equating the formation and destruction rates. These 
expressions can then be solved to give the number density 
ratios at steady state. After some algebra one obtains: 



i?(ND 3 ) 
i?(NHD 2 ) 
i?(NH 2 D) 



cR 



(4 - e) + (4 - 4e)i? 

2eR 



(4 - 2e) + (4 - 3e)R - 3ei?(ND 3 ) 
3eR 



(4 - 3e) + (4 - 2e)R - 2ei?(NHD 2 ) 



(10) 

(11) 
{12) 



Note that, although the absolute abundances will depend 
on the formation rate of NH3 , the relative abundances of 
the deuterated forms depend only on R and e, and so are 
not affected by any uncertainty regarding the kinetics of 
ammonia synthesis. 

Solutions of equations (|l0|)~(|l2|) are plotted in Figure 2. 
For R <; 1 the predicted values are proportional to R. At 
larger values of R the curves asymptotically approach the 
value determined by the branching ratios of the deuter- 
ated ammonium ions. The shaded region in Fig. 2 cov- 
ers the observed range in R; from a lower bound of 0.05 
inferred from DCO + observations (see §1), up to a max- 
imum of 0.35 implied by the N 2 D + observations of Tine 
et al. (2000). The filled circles in Fig. 2 show the observed 
NH 2 D and NHD 2 ratios; the fact that both the observed 



values imply the same underlying value of R is convinc- 
ing evidence that gas-phase chemistry is the source of the 
NH 2 D and NHD 2 in L134N. 

Note that there also exists another gas phase route to 
singly-deuterated ammonia: the reaction of N + with HD, 
analogous to reaction (§). In this case, the channel leading 
to ND + + H is favored over NH + + D, since the endother- 
micity of the former channel is only 16 K (Marquette et al. 
1988). Hence, if the H 2 ortho-to-para ratio is low so that 
reaction (^|) occurs slowly, this can lead to large enhance- 
ments of NH 2 D via the deuterated version of reaction se- 
quence (||) (Tine et al. 2000). However, this mechanism 
cannot explain the presence of doubly-deuterated ammo- 
nia. As we have shown that deuteron transfer reactions are 
more than capable of producing the observed NH 2 D and 
NHD 2 abundances, we believe that the scheme illustrated 
in Fig. 1 is sufficient to explain the ammonia deuteration 
in L134N. 



5. THE CHEMICAL STATE OF L134N 

We have shown that the observed abundances of ammo- 
nia isotopomers in L134N are consistent with the steady- 
state solutions of a relatively simple ion-molecule reaction 
scheme (see Fig. 1), if the value of R ~ 0.1. In order to 
assess the validity of our analysis, we would like to be able 
to reconcile our scheme with the observed abundances of 
simple molecular ions and their deuterated counterparts. 

The observed range of values for R implies a value of 
S of 5000-35000 (eqn. @]). In order to calculate the ex- 
pected value of S in L134N from equation (||), we need 
to know the electron fraction, x e , and 53 - kjx(mj), the 
total removal rate of H 2 D + through reactions with heavy 
species. We can calculate these values approximately from 
the the steady-state H^ concentration, which is given by 



n(H+) 



C 



Eikjx( m j) + a e (H+)x e + k f R(ED) 



and the charge conservation equation: 
x e ks x(HCO+) + .t(N 2 H+) 



x(H 3 



(13) 



(14) 
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Fig. 2. — Steady-state fractionation ratios for deuterated ammonia isotopomers, as a function of the mean deuterium 
fractionation in molecular ions, R. The curves correspond to solutions of equations (|n)|)-(]l2|) for e = 0.69, the filled circles 
show the observed values of i?(NH 2 D) and i?(NHD2), and the shaded region illustrates the range of values of R implied 
by different observations. 



Note that the rate coefficients in equation (|13|) are the 
same as in equation (||), since we assume that H2D + re- 
acts at the same rate as H^", except for electron recom- 
bination, where we use the value for a e (H^) of 1.15 x 
1CT 7 (T/300K)- - 65 measured by Sundstrom et al. (1994). 

HCO + and N2H 4 " have been observed in L134N, and 
their respective abundances are 1.2 x 10~ 8 and 7 x 10~ 10 
(Swade 1989; Womack et al. 1992a; Dickens et al. 2000). 
Adopting physical parameters of T = 10 K and n(H.2) = 
2 x 10 4 cm -3 (Dickens et al. 2000), we thus have three 
equations (0], fn§, and [Q) in three unknowns (x e , 
x(H^), and kjx(rrij)). Solving these equations for 
different values of S allows us to constrain the chemi- 
cal state of L134N. Also, because the rate coefficients 
for reactions of H^ and H2D + with heavy molecules arc 
k 2 x 10~ 9 cm 3 s _1 , we can calculate the total abundance 
of heavy molecules from our derived value of ■ kjx(mj). 

For S = 5000, we calculate x e = 1.4 x 10~ 8 , x(H+) = 
1.6 x 10 -9 , and J2j x ( m j) — 1-6 x 10~ 4 , consistent with 
little depletion of CO, N2, and O from the gas phase. On 
the other hand, for <S = 35000, we derive x e = 1.9 x 10~ 8 , 
x(H+) = 6 x 10~ 9 , and J2j x ( m j) = !- 9 x 10 ~ 5 , imply- 
ing significant depletion. For the intermediate value of 
S ~ 10 4 inferred from the deuterated ammonia fraction- 
ation, we find that partial depletion (~ 50%) of heavy 
molecules is required to account for the observed deu- 
terium enhancements. A similar conclusion has been 
reached by Roberts & Millar (2000). Finally, as our calcu- 
lated ionization levels lead to good agreement with the ob- 
served NH 3 abundance (§||), we conclude that the observed 
abundances and fractionations in L134N are well-matched 
by steady-state ion-molecule chemistry. 

6. GAS PHASE VS. GRAIN SURFACE AMMONIA 
FRACTIONATION 

Despite reservations concerning the physical conditions 
in L134N, some as yet unidentified mechanism may be 



returning mantle-formed molecules to the gas there. The 
mantle abundances of ammonia isotopomers computed nu- 
merically by Brown & Millar (1989) (see their table 2) 
effectively rule out a grain surface origin when scaled to 
the ammonia abundance in L134N. However, their sur- 
face reaction scheme permits larger D/H ratios than those 
presented, since the fractionation is proportional to the 
gas phase atomic D/H ratio, i?(D), which may be higher 
than the value assumed in their calculations. Following 
the scheme of Brown & Millar, we can derive values for 
the surface fractionation ratios: 

i? s (NH 2 D) = ^JZ(D) (15) 
i? s (NHD 2 ) = ^R(P) (16) 
i? s (ND 3 ) = ^iZ(D) (17) 

Thus, if -R(D) w 0.05, then it is also possible to explain the 
observed abundances by surface formation of ammonia. 

It may be possible, however, to discriminate between 
alternative formation mechanisms by examining the rela- 
tive scaling of the fractionation ratios. For example, grain 
surface formation implies that 

i? s (NH 2 D) : i? s (NHD 2 ) : R S (NB 3 ) = 1 : 0.33 : 0.11 (18) 

whereas, to first order, gas phase chemistry implies (using 
equations @-|||) 

i?(NH 2 D) : i?(NHD 2 ) : i?(ND 3 ) = 1 : 0.49 : 0.19 (19) 

The fact that the calculated values of R for the three iso- 
topomers always have the same ratio irrespective of their 
absolute values is apparent from Fig. 2, where the verti- 
cal separation of the curves for the three isotopomers is 
constant. 

For doubly-deuterated ammonia, the predicted ratios 
are similar, with the observed value of i?(NH 2 D) = 0.1 
implying respective values of i?(NHD 2 ) = 0.05 and 0.03 
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for gas phase and surface chemistry. Although the ob- 
served value of 0.05 agrees with our gas phase scheme, 
the observational uncertainties are too large to rule out 
surface formation. However, for triply-deuterated ammo- 
nia, the gas phase scheme predicts a fractionation almost 
twice as large; when scaled to the NH3 abundance this 
implies a ND3/NH3 ratio 2.5 times greater than the value 
predicted by surface chemistry (where we have used the 
relation ND3/NH3 = i?(ND 3 ) x i?(NHD 2 ) x i?(NH 2 D); 
cf. our definition of R in §1). With the observed val- 
ues x(NH 3 ) = 10~ 7 and i?(NH 2 D) = 0.1, we predict 
x(ND3) w 10~ n . Hence, if ND3 can be detected (or an 
upper limit determined) in L134N, it may be possible to 
determine whether the ammonia is formed in the gas or 
on the grains. 

It is worth stressing that the kind of scaling relations 
for multiply-deuterated fractionation ratios expressed by 
equations (|l^) and ([l9]) are applicable in general to all 
molecules. Thus, whereas the fractionation of singly- 
deuterated molecules reflects both the formation mecha- 
nism of the molecule and the underlying D/H ratio in the 
precursors, the relative fractionation ratios of multiply- 
deuterated molecules reflect only the formation mecha- 
nism. This fact was first appreciated by Turner (1990), 
who used the D2CO:HDCO:H2CO ratios to show that 
formaldehyde in the Orion Compact Ridge should have 
a grain surface origin. 

7. DISCUSSION 

We have shown that large abundances of NH 2 D and 
NHD 2 can be produced by gas phase chemistry in the in- 
teriors of cold dense clouds. Ammonia is deuterated via 
deuteron transfer from species such as H2D" 1 ", DCO + , and 
N2D + , followed by dissociative recombination. This mech- 
anism is able to match the observed fractionationjatios of 
both species if the underlying XD + /XH + ratio, R, equals 
0.1. 

Grain surface formation of ammonia produces distinct 
fractionation ratios, however the uncertainties in the ob- 
served abundances mean that we cannot definitively con- 
clude that deuterated NH3 is being formed in the gas. Be- 
cause the scaling of the fractionation ratios expected from 
these two processes is not the same, this raises the possi- 
bility that the NH 2 D:NHD 2 :ND 3 ratios may ultimately be 
used to determine the origin of these molecules. In particu- 
lar, this could be resolved with the detection of ND3 , which 
we predict to have an abundance of ~ 10~ n in L134N. 
Gas phase formation appears more feasible than surface 
chemistry since it is able to account for the observed NH3 
abundance without recourse to uncertain surface processes 
and desorption mechanisms. A further problem for the 
grain surface hypothesis is the fact that the large value of 
i?,(NH2D) in L134N requires a gas phase atomic D/H ra- 
tio of 0.05, but theoretical models predict an equilibrium 
value of only a few times 10" 3 at 10 K (Millar et al. 1989; 
Roberts & Millar 2000). 

We also find that the large observed molecular D/H ra- 
tios can only be reproduced if heavy elements are partially 
depleted onto grain surfaces. Therefore, it appears that 
the deuterium emission peak in L134N traces a small re- 
gion where significant amounts of CO, N2 and O are frozen 



onto grains. A similar region of enhanced D fractionation 
is known to exist in TMC-1 (Guelin et al. 1982), a dark 
cloud that appears to be physically similar to L134N. The 
fact that there appear to be several infrared sources lo- 
cated behind L134N (Snell 1981) may allow the molecu- 
lar depletion into ice mantles to be measured and crudely 
mapped. Alternatively, these spatial gradients may be 
due to the energy available when ions and neutrals have 
slightly different velocities. In this case, k r becomes the 
dominant term in the denominator of the expression for S 
(eqn. J3|), which is consequently reduced (Charnley 1998). 
Another explanation may be the existence of chemical 
bistability in interstellar clouds; Gerin et al. (1997) showed 
that molecular D /H ratios are typically reduced by a factor 
of ten in the high ionization phase steady-state solution, 
as opposed to the low ionization phase. 

Nevertheless, it seems more likely that depletion is the 
cause of the high deuteration in L134N, since the latter 
mechanisms act to reduce S whereas depletion causes S 
to increase. It is interesting to note that the observed 
range of R in L134N is always > 0.05; this is what one 
would expect for a 10 K cloud with no depletion, and is 
in fact the value derived from DCO + observations of a 
large number of cold clouds (Guelin et al. 1982; Butner 
et al. 1995). The fact that the fractionation in L134N is 
determined to be above this canonical cold cloud value is 
evidence for selective deuterium enhancement in this par- 
ticular region, as opposed to a reduction of the fraction- 
ation in the surrounding gas. On the other hand, Gerin 
et al. (1997) observed a value of R(DCO + ) = 0.003 in the 
high latitude cloud MCLD 123.5+24.9, ten times lower 
than the 'normal' value. In this instance, it would appear 
that some mechanism is indeed operating to suppress the 
D/H enhancement. 

If the depletion is particularly large, the value of R 
can become as large as unity (Millar et al. 2000). In 
this case, chemical models predict large abundances of 
HDCO and D 2 CO (Roberts & Millar 2000), so it may be 
fruitful to search for D2CO in dark clouds. Note, how- 
ever, that deuterated formaldehyde may not in fact be 
as abundant as predicted by these models where it is as- 
sumed that H2CO is deuterated via the same mechanism 
as NH3, namely deuteron transfer from H2D + , DCO + etc. 
followed by dissociative recombination. In fact, the lowest 
energy isomer of protonated formaldehyde has the form 
H2COH+, and this ion has been detected in the interstel- 
lar medium (Ohishi et al. 1996). We may therefore expect 
that deuteron transfer to H2CO will most likely result in 
H2COD + and will not lead to HDCO after dissociative 
recombination (Sen, Anicich, & Federman 1992). 

Finally, we note that the DCO + fractionation of 0.05 
observed in a number of sources by Butner et al. (1995) is 
only a factor of two less than the value of R required to 
explain the ammonia deuteration in L134N. Therefore, 
multi-deuterated ammonia may be widespread in dark 
clouds. 

Theoretical astrochemistry at NASA Ames is supported 
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